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I. INTRODUCTION 

More than twenty years ago, Berry first observed that quan- 
tum systems may retain a memory of their motion in Hilbert 
space through the acquisition of geometric phases |1]. Re- 
markably, these phase factors depend only on the geometry of 
the path traversed by the system during its evolution. Soon 
after Berry's discovery, geometric phases became a subject 
of intense theoretical and experimental studies |2]. In recent 
years, renewed interest has arisen in the study of geomet- 
ric phases in connection with quantum information process- 
ing 1 3, 4]. Indeed, geometric, or holonomic quantum com- 
putation (QC) may be useful in achieving fault tolerance |5], 
since the geometric character of thephase provides protection 
against certain classes of errors f0,'r,'K]. However, a compre- 
hensive investigation in this direction requires a generaliza- 
tion of the concept of geometric phases to the domain of open 
quantum systems, i.e., quantum systems which may decohere 
due to their interaction with an external environment. 

Geometric phases in open systems, and more recently their 
applications in holonomic QC, have been considered in a 
number of works, since the late 1980's. The first approach to 
the subject used the Schrodinger equation with non-Hermitian 
Hamiltonians 19, 10]. This is a phenomenological, non- 
rigorous approach (e.g., it cannot guarantee completely pos- 
itivity). A consistent non-Hermitian Hamiltonian description 
of an open system in general requires the theory of stochas- 
tic Schrodinger equations fllll . Nevertheless, this approach 
for the first time indicated that complex Abelian geometric 
phases should appear for systems undergoing cyclic evolu- 
tion. In Refs. flZ flsl fl4. 15, 16, 17], geometric phases ac- 
quired by the density operator were analyzed for various ex- 
plicit models within a master equation approach, but no gen- 
eral theory was formulated for open system geometric phases. 
In Refs. [8, . .18.1 , the quantum jumps method was employed to 
provide a definition of geometric phases in Markovian open 
systems (related difficulties with stochastic unravellings have 
been pointed out in Ref. ] 19]). In another approach the den- 
sity operator, expressed in its eigenbasis, was lifted to a pu- 
rified state I20l 12111 . In Ref. f22], a formalism in terms of 



mean values of distributions was presented. An interferomet- 
ric approach for evaluating geometric phases for mixed states 
evolving unitarily was introduced in Ref. ] 23] and extended to 
non-unitary evolution in Refs. L24. .25.1 . This interferometric 
approach can also be considered from a purification point of 
view ]2?, "Js"]. This multitide of different proposals have re- 
vealed various interesting facets of the problem. Nevertheless, 
the concept of adiabatic geometric phases in open systems re- 
mains unresolved in general, since most of the previous treat- 
ments did not employ an adiabatic approximation genuinely 
developed for open systems. Note that the applicability of the 
closed systems adiabatic approximation ]27] to open systems 
problems is not a priori clear and should be justified on a case- 
by-case basis. Moreover, almost all of the previous works 
on open systems geometric phases were concerned with the 
Abelian ^Be i ry p hase) case. Exceptions are the very recent 
Refs. ]8, 28,^9], which discuss both non-adiabatic and adi- 
abatic dynamics, but employ the standard adiabatic theorem 
for closed systems in the latter case. 

In this work, we introduce a self-consistent open-systems 
framework, based on a recent generalization of the adiabatic 
approximation fsoil . which allows for a general definition and 
evaluation of both Abelian and non-Abelian geometric phases 
in open systems undergoing cyclic adiabatic evolution. As 
we shall show, this approach yields new insights and lends it- 
self to a simple and elegant generalization of the concept of 
geometric phases. An important feature emerging from this 
picture is the appearance of a distinguished time-scale for the 
observation of adiabatic geometric phases in open systems. 
We illustrate our results by considering the canonical exam- 
ple of a spin- 1/2 in a magnetic field. In this example, we 
find a remarkable robustness of the geometric phase against 
both dephasing and spontaneous emission in the instantaneous 
renormalized Hamiltonian eigenstate basis. This result should 
have a positive impact on the robustness against external dis- 
turbances of holonomic QC. 

We note that an alternative theory for adiabaticity in open 
systems was recently developed by Thunstrom, Aberg and 
Sjoqvist, for systems coupled weakly to their environment. 
This theory was then employed to study a non-Abelian geo- 
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metric phase gate in holonomic QC f2^. The main differ- 
ence between the approach of Ref. |26] and our approach to 
adiabaticity in open systems, is that Thunstrom et al. empha- 
size the decoupling of eigenspaces of the system Hamiltonian, 
while we focus on the entire superoperator. Central conclu- 
sions, such as the breakdown of adiabaticity in open systems, 
are shared by the two approaches. A full comparison is be- 
yond the scope of the present work. 

The structure of this paper is as follows. In the next section 
we review the coherence vector approach to solving quantum 
master equations, the Jordan form of the superoperator, and 
the notion of adiabaticity in open quantum systems. In Sec- 
tion Hn] we derive the geometric phase in open systems, in 
both the Abelian (Berry phase) and non-Abelian cases. We 
prove that this geometric phase satisfies the expected prop- 
erties, such as a proper closed system limit and gauge in- 
variance. In Section Hvl we focus on applications, namely, 
we show that our theory predicts that there is a distinguished 
time-scale for open system geometric phases (related to the 
breakdown of adiabaticity), and then consider the example of 
a spin- 1/2 coupled to a slowly varying magnetic field, in the 
presence of dephasing and spontaneous emission. We con- 
clude in Section El 

II. MASTER EQUATIONS AND THE ADIABATIC 
REGIME OF OPEN QUANTUM SYSTEMS 

Open quantum systems typically do not undergo unitary dy- 
namics, i.e., they are not governed by the Schrodinger equa- 
tion, or even by its non-Hermitian generalization I^ Uoll . In- 
stead, quite generally we may consider open quantum systems 
evolving under a convolutionless master equation 13 1.1 

dp/dt = C[R{t)]p{t), (1) 

where £ is a superoperator which depends on time only 
through a set of parameters R{t) = R. The Lindblad equa- 
tion 13211 is an important example of this class of master equa- 
tions: 

dp/dt = -I [H, P] + \J2 (F- P^l^ + F*/^' rll) . (2) 

i 

where we have suppressed the explicit dependence of the 
operators on R{t). Here H is the effective Hamiltonian of 
the open system (it is renormalized, i.e., contains the "Lamb 
shift" - the unitary contribution of the system-bath interac- 
tion the Ti are operators describing the system- 
bath interaction, and we work in h = 1 units. In this work 
we consider the general class of convolutionless master equa- 
tions 0, and in a later section illustrate our formalism with 
an example using the case of Eq. In this example, of a 
spin-1/2 in a magnetic field, we allow both H and the Ti to 
depend on R{t). In a slight abuse of nomenclature, we will 
refer to the implicitly time-dependent generator C [Eq. Q] as 
the Lindblad superoperator and the Vi [Eq. (|2j] as Lindblad 
operators. This terminology is usually associated with time- 
independent generators 132. .331,1 . but recent work has clarified 



the conditions under which Eq. (|2} with time-dependent 
can be derived in the usual Davies L35i] weak-coupling limit 
Js^. What is important to note is that the microscopic weak- 
coupling limit derivation is consistent with the general class 
of master equations postulated here, wherein the Lindblad op- 
erators depend implicitly on time through their explicit de- 
pendence on external control fields. More specifically, in the 
microscopic derivation one shows that the Lindblad opera- 
tors are the Fourier components of the time-dependent sys- 
tem operator terms in the system-bath interaction Hamilto- 
nian, where the time-dependence arises by working in the 
interaction picture with respect to the renormalized system 
Hamiltonian |36]. This is how the F^'s appearing here and 
in our example below must be interpreted. 

In the superoperator formalism, the density matrix for 
a quantum state in a D-dimensional Hilbert space is rep- 
resented by a D^-dimensional "coherence vector" \p)) = 
{pi, P2, ■ ■ ■ , Pd'^) and the Lindblad superoperator £ becomes 
a X Z?'^ -dimensional supermatrix |33], so that the master 
equation Q can be written as linear vector equation in D^- 
dimensional Hilbert-Schmidt space, in the form 

d\p))/dt = C[Rmp)). (3) 

Such a representation can be generated, e.g., by introducing a 
basis of Hermitian, trace-orthogonal, and traceless operators 
[e.g., the ZJ-dimensional irreducible representation of the gen- 
erators of su(I?)], whence the pi are the expansion coefficients 
of p in this basis |33], with pi the coefficient of / (the identity 
matrix). In this case, the condition Trp^ < 1 corresponds to 
|||p))|| < 1. P = to Pi = p*, and positive semidefiniteness 
of p is expressed in terms of inequalities satisfied by certain 
Casimir invariants [e.g., of su(£')] 1371 l38(l . A simple and 
well-known example of this procedure is the representation of 
the density operator of a two-level system (qubit) on the Bloch 
sphere, via p = {I2 + ^ ■ ^)/2, where if ~ {a^, (Ty,az) is 
the vector of Pauli matrices [generators of su(2)] and I2 is the 
2x2 identity matrix. 

The master equation generates a non-unitary evolution 
since C is non-Hermitian. In fact, C need not even be a nor- 
mal operator (C^ C ^ CC^). Therefore C is generally not di- 
agonalizable, i.e., it does not possess a complete set of lin- 
early independent eigenvectors. Equivalently, it cannot be 
put into diagonal form via a similarity transformation. How- 
ever, one can always apply a similarity transformation to 
C which puts it into the (block-diagonal) Jordan canonical 
form 1 39], namely, Cj = S^^CS. The Jordan form Cj 
of a X matrix £ is a direct sum of blocks of the 
form £j = ®^=iJa (Q! enumerates Jordan blocks), where 
771 < is the number of linearly independent eigenvectors 
of £, X]^=i "-a — where = dim Ja is the dimension 
of the ath Jordan block, and = Xaln^ + Ka where Aq. 
is the ath (generally complex-valued) Lindblad-Jordan (LJ) 
eigenvalue of C (obtained as roots of the characteristic poly- 
nomial), In^ is the Ua X Ua dimensional identity matrix, and 
Ka is a nilpotent matrix with elements {Kajij — <5i.j-i (I's 
above the main diagonal), where S is the Kronecker symbol. 
Instantaneous right {\vj^\R{t)]))} and left {{{£i'\R{t)]\} 
bi-orthonormal bases in Hilbert-Schmidt space can always 
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be systematically constructed such that they obey the or- 
thonormahty condition {{£^a\vf)) = S^pS*^ fs^. Here 
superscripts enumerate basis states inside a given Jordan 
block e {0,...,na — 1}). When C is diagonalizable, 
{\v'-f [R{t)]))} and {((£L'\R{t)]\} are simply the bases of 
right and left eigenvectors of £, respectively. If C is not di- 
agonalizable, these right and left bases can be constructed by 
suitably completing the set of right and left eigenvectors of £ 
(See Ref. |30| and also Appendix lAlfor a detailed discussion 
of the left and right basis vectors, including their complete- 
ness relation). Based on the above considerations we gave, 
in Ref. ll30ll . a definition of adiabaticity in open quantum sys- 
tems: 

An open quantum system is said to undergo 
adiabatic dynamics when its Hilbert-Schmidt 
space can be decomposed into decoupled LJ- 
eigenspaces with distinct, time-continuous, and 
non-crossing instantaneous eigenvalues of C. 

Note that the key to establishing the concept of adiabaticity in 
open systems is to replace the idea of decoupling of the Hamil- 
tonian eigenstates by that of decoupling of the Jordan blocks 
of the Lindblad superoperator (each Jordan block is associ- 
ated with an independent eigenstate of C). The definition of 
adiabaticity given above implies a condition on the total evolu- 
tion time T which generalizes the well-known closed-systems 
condition, T > maxo<s<i |(fc(s)|^^^^|m(s))|/|5„fe(s)p, 
where s = t/T is the normalized time, H is the time- 
dependent Hamiltonian, |fc) and |m) are eigenstates of H, 
and gmk is the energy gap between these two states. For fur- 
ther details we refer the reader to Refs. the con- 
dition in the case of 1 -dimensional Jordan blocks is given in 
Eq. i32\ below. The theory developed in Ref. 1 30] applies also 
in the more general case of explicitly time-dependent gener- 
ators C{t), but since we focus here on geometric phases, we 
shall only consider imphcit time-dependence as in C[R{t)]. 



Ul. GEOMETRIC PHASES FOR OPEN SYSTEMS IN 
CYCLIC ADIABATIC EVOLUTION 

In order to define geometric phases in open systems, we 
expand the coherence vector in the instantaneous right vector 

basis {\'Df[R{t)]))}as 

13=1 j=0 

where we have explicitly factored out the dynamical phase 
exp[J'p Xi3{t')dt']. The coefficients {p^p\t)} play the role of 
"geometric" (non-dynamical) amplitudes. We assume that the 
open system is in the adiabatic regime, i.e., Jordan blocks as- 
sociated to distinct eigenvalues evolve in a decoupled manner 
(recall the definition of open sytems adiabaticity given above). 



Then, Eqs. (|3},(E} together yield 

E Y.^^^^"^))- (5) 

A condition on the total evolution time, which allows for the 
neglect of coupling between Jordan blocks used in deriving 
Eq. was given in Ref. |30]. It is reproduced for the case 
of 1 -dimensional Jordan blocks in Eq. (I32> below. Note that, 
due to the restriction \p = \a, the dynamical phase has disap- 
peared. For closed systems, Abelian geometric phases are as- 
sociated with non-degenerate levels of the Hamiltonian, while 
non-Abelian phases appear in the case of degeneracy. In the 
latter case, a subspace of the Hilbert space acquires a geo- 
metric phase which is given by a matrix rather than a scalar. 
Here, for open systems, one-dimensional Jordan blocks can 
be associated either with Abelian or non-Abelian geometric 
phases (depending on the possibility of degeneracy) while 
multi-dimensional Jordan blocks are naturally tied to a non- 
AbeUan phase. 

A. The Abelian Case: Generalized Berry Phase 

Consider the simple case of a non-degenerate one- 
dimensional Jordan block (a block that that is a 1 x 1 subma- 
trix containing an eigenvalue of C). In this case, the absence 
of degeneracy implies in Eq. (|5} that ~ \a ^ oi = (3 
(non-degenerate blocks). Moreover, since the blocks are as- 
sumed to be one-dimensional we have tIq, — 1, which allows 
us to remove the upper indices in Eq. (|5}, resulting in pa — 
—Pa{{£a\T^a))- The solution of this equation is Pa{t) — 
p„(0)exp [i7„(t)], with7„(i) = i J*{{£a,{t')\t)a,{t')))dt'. In 
order to establish the geometric character of (t) we now re- 
call that C depends on time implicitly through the parameters 
R{t). Then, for a cyclic evolution in parameter space along a 
closed curve C, we obtain that the Abelian geometric phase 
associated with the Jordan block a is given by 

7„(C) =i£((£„(i?)|vPa(i?)» -rfi?. (6) 

This elegant generalized expression for the geometric phase, 
which bears similarity to the original Berry formula |[l|l, is 
our first main result. As expected for open systems, 7^ (C) is 
complex, since {{£a\ and IPq)) are not related by transpose 
conjugation. Thus, the geometric phase may have real and 
imaginary contributions, the latter affecting the visibility of 
the phase. 

In Refs. I?, 1^ Garrison and Wright, and Dattoli et al., 
found an expression for the open-systems Berry phase that 
resembles our Eq. (|6j. Their result is 

7„(C) = I £{e^{R)\^\MR)) ■ dR. (7) 

Here and are a bi-orthonormal set of eigenvec- 

tors of a non-Hermitian Hamiltonian H and its Hermitian con- 
jugate iJ^ respectively. There are some important method- 
ological and technical differences between this and our result. 
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First, here, instead of working with a phenomenological non- 
Hermitian Hamiltonian, we started from the outset with a fully 
consistent master equation approach, where the left {((fall 
and right {II'q))} basis vectors are associated with the dy- 
namical superoperator £, rather than with the non-Hermitian 
Hamiltonian. Second, as a result in our case, the basis vec- 
tors span the -dimensional Hilbert-Schmidt space, whereas 
in non-Hermitian Hamiltonian case the geometric phase ex- 
pression involves vectors in the usual _D-dimensional Hilbert 
space. As we note below, this implies that in our case is 
a relative, not absolute, geometric phase, and hence there is 
in general no connection between the expressions (|6|l and Q. 
Third, unlike Refs. ^ JO.] , where adiabaticity is imported 
from the theory of closed systems, we work within a consis- 
tent theory of adiabaticity for open systems, as formulated in 
Ref. \M- 

The expression for 7a (C) exhibits a number of important 
properties expected from a good definition of a geometric 
phase: 

• Geometric character. ^a{C) is geometric, i.e., it de- 
pends only on the path traversed in parameter space. 



pect in which our expressions differ from the ones de- 
rived using non-Hermitian Hamiltonians i^ fToll . 



B. The non-Abelian Case: Generalized Holonomic Connection 

Let us now generalize these considerations to degenerate 
one-dimensional Jordan blocks, whence the geometric phase 
becomes non-Abelian. From Eq. ^ we obtain 



/3| A^=A„ 



(8) 



Each decoupled subspace is associated with a different value 
of Aq, and is spanned by the set {\'Di3{R))) | A^j ~ A^}. 
Then, enumerating this list for each decoupled subspace (de- 
noted by Aq) as {V^ 



(1) 



)),\V^^;>)),...,\VZ')) (the left basis 



{{{£a{R)\} is similarly enumerated), with G the degener- 
acy (dimension of the decoupled subspace), we have that 



Ef=iPA ^(4'^I^A^»- Writing this equation in 



a vector notation, we obtain 



• Gauge invariance: ^a{C) is gauge invariant, i.e., we 
cannot modify (or eliminate) the geometric phase by 
redefining {{£a\ or |I?q)) by multiplying one of them 
by a complex factor x{R) exp [iv{R)]. Indeed, let us 
define = xexp (iz.)|I?a)) (x(i?) ^ Vi?). Re- 
definition of right-vectors automatically implies redefi- 
nition of left- vectors due to the normalization constraint 
{{£a\'Da)) = 1, so that ((541 = ((f„|x-iexp(-«z/). 
Therefore, ((f^lVPL)) - ((^^Ivpa)) + (Vx)/x + 
iS/v. Gauge invariance then follows from the computa- 
tion of 7^ using Eq. (|6j, with Stokes's theorem leading 
to 7q(C) — ^ci{C). Below we provide a detailed proof 
of gauge invariance in the non-Abelian case, which in- 
cludes the latter as a special case. 

• Closed system limit: if the interaction with the bath 
vanishes, 7q(C) reduces to the usual difference of 
geometric phases acquired by the density operator in 
the closed case. In order to prove this, consider the 
expansion of the vectors and {{£a\ in a basis 
{lD,J^i\i = 1,...,!?^ - 1}, where Id is the x 
identity matrix and are traceless Hermitian matri- 



ces, with Tr(AiAj 



Then, by using the normal- 



ization condition {{EoXDa)) = 1 and the matrix inner 
product {{u\v)) = [1/ D)Tx{u^v) js^, we obtain in 
the closed-case limit |I?q)) \/15\ipm){'^n\ [equiva- 
lently, {{Ea\ -> \fD\i'n) {i^mW, with {|V'm>} denoting 
a set of normalized eigenstates of the Hamiltonian op- 
erator. Therefore, Eq. (|5Jl yields 7^ ^f^"^ — 
which is exactly the difference of phases acquired by 
the density matrix in closed systems. Note that only 
phase differences are experimentally observable, so that 
the fact that our expression for the geometric phase in- 
volves phase differences, rather than an absolute phase, 
is natural. As mentioned above, this is an important as- 



-(1a„ •^)Pa 



(9) 



where P^^ — (^P^\\ ■ ■ ■ , P^x'^^ is a vector in Hilbert- 
Schmidt space (superscript t denotes transposition) and 



lvl<^» 



((4'jlv|2?i!^)> \ 



V((4!'ivi<')) ••• ((4!'ivPl!^» ; 



7(G) I 



(10) 

Note that each element of yl is a vector in parameter space 
(we use bold-face and arrow superscripts to denote vectors 
in Hilbert-Schmidt space and parameter space, respectively). 
The non-Abelian geometric phase in a cyclic evolution asso- 
ciated with a degenerate level Aq, is determined by the solu- 
tion of Eq. (|9j, which is formally provided by Pa„ (C) = 
U Pa„ (0), where 



U^Ve 



- Ax 



■dR 



(11) 



is the corresponding Wilson loop, and V denotes path- 
ordering. 

Equations ilO\ and (II 1> constitute our second main result. 
They are the generalization of the concept of non-Abelian ge- 
ometric phases to the open systems case. In particular, the 
matrix [Eq. ilQ \] naturally generalizes the Wilczek-Zee 
gauge potential [40], also known as the holonomic connec- 
tion. 

A non-Abelian geometric phase will also appear in the case 
of multi-dimensional Jordan blocks. However, in this case, it 
is not possible to obtain a general analytical solution due to the 
presence of the term pi*^^^ in Eq. (|5}- One should then solve 
Eq. (|5} on a case by case basis for all pairs (a, i). This yields 
a set of coupled differential equations in a ladder structure. 
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The geometric character of the non-Abelian geometric 
phase is evident from the expression il 1> for the Wilson 
loop operator (it depends only on the path and not on its 
parametrization). The closed system limit is obtained in a 
manner exactly analogous to the proof above for the Abelian 
case. What is left, therefore, in order to demonstrate that our 
expressions for the generalized non-Abelian geometric phase 
have the desired properties, is a proof of gauge invariance of 
the eigenvalues of the Wilson loop (this is the same as in the 
closed systems case |40|, where the Wilson loop itself, with- 
out taking the trace out, is not gauge invariant). We consider 
this issue next. 



C. Gauge invariance of the non-Abelian geometric pliase in 
open systems 

Let us show that the geometric phase derived in Eq. (^Ji is 
invariant under gauge transformations. First, we can rewrite 
Eq. (|9jl in the form 

i^--f:^i'pt (12) 

where 

4^) ^ = T.(4%^ = (4?) • ^, 

^ (13) 

where we have introduced the explicit dependence on R{t) — 
{xi{t), ...,Xn{t)), and 

(14) 



Let us now apply a local gauge transformation n{R{t)) on the 
right eigenvectors, 

d 

where fl is an arbitrary complex matrix. Here is taken as a 
function of R{t) instead of t because the basis vectors depend 
on t only through the parameters R{t). The normalization 
constraint {{£x^\T^'\^)) = 1 implies 

((4!'i = E((4:^i"cV- (16) 

c 

The gauge transformation of A^^"''^ = ((f^''-* |I>^''-')), is then 
given by 

cd 

= {nA{^n-^)ab + inn-\b, (i7) 

where superscript t denotes transposition. Therefore 

A'l^ = fiAl^n^^ + hn^^, (18) 

which proves that transforms as a gauge potential. 

Now let us show that the Wilson loop has gauge-invariant 
eigenvalues (a similar proof for the closed case can be found 
in Ref. |41]). As a first step, consider the Wilson operator for 
an open curve 



U = Pexpl-y^ A-dRj ^Vexpi^~J dt' A^{x) — 



I- I dtiA^,{ti)-^ + dtiAf,,{h)-^ dhA^,{t2)^ + ... , (19) 



where repeated indices are summed over, and we suppress the subscript for notational simplicity. The transposed Wilson 
operator then yields 



Note the inversion of the order of the operators due to the gauge invariance of Eq. (HH . After a gauge transformation, 
transposition. Therefore, the transposed Wilson operator Eq. (I21t reads: 
W = U* obeys the differential equation 

dW dr^^ 

^^ + WA-/^^0^ (21, ^ + »"<^-»- <^^' 



We can determine the gauge transformation of W by imposing 



where primes indicate gauge-transformed operators. Note that 
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A^, which is given by 

= (23) 

transforms, according to Eq. ( I18> . under a gauge transforma- 
tion as 

A'^ ^ nAl^-^ + {df,n)n-\ (24) 

Then, using Eq. (|23, we obtain 

Since il is arbitrary it follows from Eq. j25t that gauge invari- 



ance of the equation of motion implies 

w^w'^nwn^\ (26) 

where fl^^ = ^l^^{x^{t)) and fl is independent of x^{t). 

The gauge transformation of a product of paths allows us 
to further restrict (l: we can show that il = f2~^(xo), where 
here and below Xi = a;''(ti) [x^{to) is the initial position 
R{0))]. To see this, consider an open curve F = Fi + F2, 
where Fi is a continuous curve in the interval [xq, Xa] and F2 
is a continuous curve in the interval [xq, xt]. The transposed 
Wilson operators Ty(Fi) and M^(F2) associated with these 
curves are, according to Eq. i20\ . given by 



dxt"^ 



"^1 

* (f.\ 

dti 



dhAlih)-^ + 



dti\ dhAl^{t2)^Al^ih)^ + 



dti 



(27) 



Then, under gauge transformation, we have from Eq. J26l i: 



IV. APPLICATIONS 



VK(Fi) ^ niW{Ti)n{xa)-^ 
VF(F2) ^ n2WiT2Mxb)-' 
w{T) nw{r)n{xb)~'^. 



(28) 



Then, by applying a gauge transformation on VF(F) = 
VF(Fi)H/(F2), we obtain 

nw{r)n{xb)-^ = niW{Ti)n{xa)-^n2W{T2Mxb)-'. 

(29) 

This implies (i) O2 = n{xa), i.e., ft depends only on the 
initial position, and (ii) J7 = J7i ($7 must be independent of 
the path index). We can therefore write (l = r2(a;o).Thus, 
for a general open curve F in the interval [xQ,Xa], the 
transposed Wilson operator transforms as — > W' = 
il{xo)Wft~^{xa)- For a closed curve C, we then get W — > 
W' — il{xo)Wfl^^ (xq), or, expressed in terms of the Wilson 
loop itself: 



(30) 



Thus, the Wilson loop transforms as a similarity transforma- 
tion, and consequently its eigenvalues are gauge invariant, as 
desired. In particular, this implies that the gauge-transformed 
amplitude vector P';^ = (fi*)^^ (xq) Pa^ obyes the same ge- 
ometric phase transformation rule as Pa„ : 



= (fl*)-i(xo)PA„(C)=P' (C). 



(31) 



A. A Distinguished Time-Scale for Open System Geometric 
Phases 



As a first application of these general considerations we 
now show the existence of a distinguished time-scale for the 
observation of open-system geometric phases. To this end, it 
is convenient to express the variables in terms of the dimen- 
sionless time s ~ t/T, where T denotes the total evolution 
time. Then, adiabatic dynamics in the interval < s < 1 
occurs if and only if the following time condition is satis- 
fied: T ^ maxQ,{T^}, where denotes the crossover time 
for the Jordan block Jq, ll30ll . For the particular case of one- 
dimensional blocks, we have I30ll42ll 



= max 

0<s<l 



(32) 



where 



^f3a{s) = / UJfjais)ds 

Jo 

W/3a(s) = 7/3(5) - 7a (s) 

is the gap between Jordan eigenvalues, 

V,^{s)=p,{s){{£^{s)\^\Vpis))) 



(33) 



(34) 
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are the matrix elements of the time-derivative of the Lindblad 
superoperator, and 



Qf3a{s) = Vi3^is)/ujfiJi 



(35) 



Note that a quantity analogous to Qpa, namely the time- 
derivate of the Hamiltonian divided by the square of the spec- 
tral gap, appears in the standard condition for adiabaticity 
in closed systems fjflll . In the expression for V3a(s), upper 

indices in p^p [s) and in the basis vectors {jP^''-' (s)))} and 

{ ( (£q ■* (t) I } were removed because the Jordan blocks are one- 
dimensional. The crossover time provides a decoupling 
timescale for each Jordan block: provided T ^ the Jor- 
dan block J a is adiabatically decoupled from all other blocks 
associated to a dijferent eigenvalue. The general expression 
for in the case of multi-dimensional Jordan blocks, as well 
as a more detailed discussion of its meaning, are given in 
Refs. SEl- 

Now, the important observation is that the decoupling time- 
scale is finite due to the presence of complex exponentials in 
llsoll . which have real and imaginary parts in the case of 
open systems. Therefore, since the geometric phases are de- 
fined in the adiabatic regime, they will only be observable 
during the finite time in which the Jordan blocks are decou- 
pled. This fact implies the existence of a distinguished, finite 
time-scale for geometric phases in open systems. Such a time- 
scale was noted in Ref. Ill6(l in the context of a specific exam- 
ple, namely the case of a spin- 1/2 particle in a magnetic field. 
Finite adiabaticity time-scales have been revealed as a gen- 
eral property of open systems fl^ l3(l l43ll . a fact which has 
also been observed in adiabatic QC, both theoretically B2il 
and experimentally |44]. Physically, the reason for this phe- 
nomenon is the broadening of the system energy levels due to 
the presence of a dense spectrum of bath energy levels, until 
the broadened system energies overlap. When this happens 
different eigenspaces may no longer be decoupled (provided 
there are no selection rules preventing the coupling), and the 
adiabatic approximation breaks down. In the case of static 
Hamiltonians this is known as quantum diffusion 1 45 1 . 



Spin 1/2 in a time-dependent magnetic field under 
decoherence 



As an illustration of the general theory presented above, 
let us consider the canonical example of a spin-1/2 in a 
time-dependent magnetic field, originally considered by Berry 
in the context of closed quantum systems |1]. The renor- 
malized Hamiltonian of the system is given by H{B) = 
— /iS* • B, where S = (l/2)(cr^, cr^, cr^) is the spin op- 
erator, with CT* (i — X, y, z) denoting the Pauli matrices. 
Bit) = (Bxjt), By(t), Bz(t)) is a time-dependent magnetic 
field (including the Lamb shift correction 11331 0411 '). and /i is 
a constant. A standard evaluation of the geometric phase in 
this case yields 7§!°"«<^(C) = ±Vl{C)/2, where 7§!°"'"*(C) 
are the geometric phases associated with the energy levels 
E± = ±(^/2)B, with B = \B{t)\, and n{C) being the solid 



angle subtended by the closed curve C traversed by the mag- 
netic field in parameter space. 

In the weak-coupling regime, it is common to consider 
decoherence in the eigenbasis of the renormalized system 
Hamiltonian ll36l l46l Wh . Let us now analyze the effects 
of decoherence in this basis, assuming that open-systems dy- 
namics is described by the master equation We consider 
two important sources of decoherence, namely, dephasing and 
spontaneous emission in the eigenenergy basis. The Lindblad 
operators modelling these processes are given, respectively, 
by = P^W{B)<7,W^B) and r_ = |3.W{B)a^W^B) 
(ct_|1) = 2|0),cr_|0) = 0), where (3:, and (3. are the er- 
ror probabilities per unit time and W{B) is the unitary ma- 
trix which diagonalizes H{B). The Lindblad superoperator is 
then given by C{B) = H(-B) + 7^(i?), where n{B) is the 
Hamiltonian superoperator [obtained from H{B)] and Tl{B) 
is the superoperator containing the decoherence contribution 
[obtained from and F_]. In this case explicit calculation 
reveals that (i) H and TZ are diagonalizable, (ii) [H, TZ] = 0. 
Hence C, Ti and TZ have a common eigenstate basis, and in 
particular it follows that C has only one-dimensional Jordan 
blocks, whence it is diagonalizable. Thus, bearing in mind 
that the eigenstate basis for 7i(i?) is independent of (3z and 
/?_, it follows that the eigenstate basis for C{B) is also in- 
dependent of Pz and /3_ . This implies that the adiabatic ge- 
ometric phases, which can be computed here from Eq. 
are robust against dephasing (F^) and spontaneous emission 
(F_). In this case, integration in parameter space, which is the 
relevant space for adiabatic geometric phases, is not affected 
by decoherence. Nevertheless, the adiabaticity crossover time 
does depend on (3^ and /3_ through the eigenvalues of C. 

The robustness against dephasing is in agreement with 
Ref. [18], but obtained here in a totally different framework. 
A microscopic derivation of the geometric phase for a spin- 
1/2 in a magnetic field was developed in Ref. [17], with no 
robustness against dephasing detected. However, note that 
the robustness of the geometric phase obviously depends on 
the basis in which the environment acts. As is clear from 
Eq. (4) of Ref. 1 17], Whitney et al. consider dephasing in 
a fixed (time-independent) basis, where robustness is absent. 
In this basis, our approach is in agreement with this lack of 
robustness. However, in the example we discuss here, we 
consider decoherence in the instantaneous eigenenergy basis. 
In the weak-coupling regime, one should consider decoher- 
ence in the eigenbasis of the system Hamiltonian (see, e.g., 
Refs. .47.1 ). This follows from the fundamental Davies 
derivation of the quantum Markovian master equation ]35l, 
which was recently reviewed and generalized in a context 
relevant to ours in Ref. [36] , and which shows that in the 
Markovian limit time-dependent system Hamiltonians are al- 
ways coupled to the Lindblad operators. This difference in 
basis explains the apparent disparity between Ref. ]17] and 
our result. Then, in the instantaneous eigenbasis, we obtain 
robustness against both dephasing and spontaneous emission 
for adiabatic evolution as a simple consequence of the com- 
mutation relation between the Hamiltonian superoperator Ti 
and the corresponding decoherence superoperators TZ. In fact. 
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the commutation between H and TZ provides a general suffi- 
cient condition for robustness of adiabatic geometric phases 
against TZ. Note also that the non-adiabatic geometric phase 
is usually affected by corrections due to the system-bath inter- 
action, in particular in the case of spontaneous emission fT3l . 
Remarkably, by imposing adiabaticity on the open system, ro- 
bustness of the geometric phase against this decoherence pro- 
cess is obtained. 




FIG. 1: Maximum value of the ratio T^/T of the crossover time to 
the total evolution time, taken over all the Jordan blocks, as a function 
of T for a spin-1/2 particle in a magnetic field undergoing dephasing. 
Parameter values are given in the text. The adiabatic interval requires 
Ta <C T. Observe that the stronger the magnetic field (at fixed 
decoherence rate), the better the adiabatic approximation, and hence 
the longer we can observe the geometric phase. 

We stress that, to the best of our knowledge, the approach 
presented here for dealing with geometric phases is the first 
to predict robustness against both dephasing and spontaneous 
emission for a spin-1/2. We expect that such a robustness will 
serve as a useful protection mechanism in holonomic QC (see, 
e.g., Ref. [8] for difficulties in the correction of spontaneous 
emission). The robustness is reminiscent of the emergence 
of a decoherence-free subspace (DFS) i48il . but unlike the 
symmetry-driven appearance of the latter, here the robustness 
is due to a, fundamentally different, adiabatic mechanism. Re- 
lated observations were made, using very different methods, 
in Ref. |49|, for a system in a DFS, coupled to a cyclicly 
evolving reservoir 

As for closed systems, B plays an important role in set- 
ting the time interval for the observation of the geometric 
phase. The reason is that the Zeeman effect further splits the 
system energy levels, thus postponing the breakdown of adi- 
abaticity due to overlap caused by environmentally induced 
broadening. This behavior is illustrated in Fig.^for the case 
of dephasing, where we take f3z — 0.1 (in units such that 
fi = 1) and use the following spherically symmetric config- 
uration for the magnetic field: Bx{s) — i?cos(27rs) sin6', 
By{s) ~ S sin(27rs) sin 6*, and B^ = Bcosd, with 6 de- 
noting the azimuthal angle, set at 7r/3. The initial state of 
the system is chosen to be an equal superposition of energy 
eigenstates. Due to the commutation relation [HjTZ] — 0, 




-1.06 1 ' 1 ' 1 ' 1 ' 1 ' 1 

0.1 0.2 0.3 0.4 0.5 

K 

FIG. 2: Geometric phase (in units of tt) of a spin-1/2 in a magnetic 
field, as a function of the decoherence strengths /3„ (n is chosen as 
z, — , or X depending on the decoherence process, as indicated in the 
legend). The magnetic field is applied in a spherically symmetric 
configuration whose parameters are chosen as in Fig. □ The ge- 
ometric phase 7 \^{C) plotted is for the Jordan block associated 

~ in the limit of vanishing decoherence - with the closed-system 
geometric phase difference 7l'°='=d(c) _ Y^°^^<i(C), which equals 
-27r[l - cos [0 = 7r/3)] = -tt. 



we find that the four Jordan blocks are associated with the set 
of Hilbert-Schmidt space vectors {\ipm) {tpn\} {m,n = 1,2), 
where {li/"™)} are the normalized eigenstates of the system 
Hamiltonian. Then, direct computation of 7^ from Eq. (|6jl 
(for s = 1) yields ja = ±27rcos0 for the vectors |'0m)('0n| 
with m n. Up to an unimportant 27r factor, these are exactly 
the differences of geometric phases ±27r(l — cos 9) appearing 
in the density operator for these states in the closed case. As 
suggested from the analytical treatment above, similar results 
hold for spontaneous emission. 

It should be noticed, however, that robustness is, naturally, 
not universal; e.g., it does not hold for the bit-flip channel 
Tx{B) — PxW{B)axW{By , since does not commute 
with the Hamiltonian superoperator This result is illustrated 
in Fig.|2j where it is shown that the real part of the geomet- 
ric phase under bit-flip is slightly affected by the decoherence 
process. This is in contrast with the robustness of dephasing 
and spontaneous emission. The corresponding imaginary part 
of the geometric phase (not shown) has negligible variation, 
which means that the visibility of the phase is not significantly 
affected by the bit-flip channel for the decoherence strengths 
considered in the plot. 

V. CONCLUSIONS 

We have introduced a general framework for geometric 
phases in open systems undergoing cyclic adiabatic evolution. 
Expressions which naturally generalize the familiar closed 
systems' (Abelian) Berry phase and (non-Abelian) Wilczek- 
Zee gauge potential and Wilson loop were derived, and their 
gauge invariance proven. An important feature of our ap- 
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proach is the existence of a distinguished time-scale for the 
observation of the adiabatic geometric phase. This prop- 
erty imposes time constraints on reahstic schemes for holo- 
nomic QC based on adiabatic phases. Remarkably, robust- 
ness against dephasing and spontaneous emission was found 
for the geometric phase acquired by a spin- 1/2 in a magnetic 
field. This robustness is, however, not universal; e.g., it does 
not hold for the bit-flip channel, since the latter does not com- 
mute with the Hamiltonian superoperator. The development 
of methods for overcoming decoherence affecting geometric 
phases is therefore of significant interest L 8. .50.1 . 
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APPENDIX A 

We define the right { |X>^-''^ [R{t)])) } and left {{{Si'^ | } 

basis vectors associated with £[R{t)] and prove their bi- 
orthogonality and completeness. Instantaneous right and left 
eigenstates of a general time-dependent superoperator C{t) 
are defined by 

-CWI^aW)) = Xait)\ram, (M) 

((Q„(i)|£(t) = ((Q„(t)|A„(0, (A2) 

where possible degeneracies correspond to Aq, = Xp, with 
a f3. In other words, we reserve a different index a for each 
independent eigenvector since each eigenvector is in a distinct 
Jordan block. It follows from Eqs. (IA1> and (IA2> that, on the 
one hand 

{{Qpimm \Vc.m) = KmiQ^it) ip.w)), (A3) 

while on the other hand this equals 

mmm)\Vo.m ^ {{Q0{t)\v^m Mt)- (A4) 

Therefore for 7^ A^g, we have {{Qait)\Vp{t))) = 0. 

The left and right eigenstates can be easily identified when 
the Lindblad superoperator is in the Jordan form Cj{t) — 
S-\t)CSit). Denoting iVc^m.; - S-\t) \Vc.m, i.e., 
the right "Jordan basis" (note the J subscript) eigenstate of 
Cj{t) associated to a Jordan block J^, then Eq. (IA1> implies 
that \Va{t)))j is time-independent and, after normalization, 
is given by 



(1,0,0, 



,0 



(A5) 



where only the vector components associated to the Jordan 
block Ja are shown, with all the others vanishing. In order to 
have a complete basis we shall define new states, which will 



be chosen so that they preserve the block structure of Cj{t). 
A suitable set of additional vectors is 



= 0,1,0, ...,0 



I = (0,,...,0,l)*, 



(A6) 



where again all the components outside Ja are zero. This sim- 
ple vector structure allows for the derivation of the expression 

Cjit) \vi^^))j = \vl^-'^))j + Kit) \vi^^))j, (A7) 



with = \Va))j and \V'a "')).! = 0. The set 

||I?Q j = 0, {ua ~ 1)} immediately be related 
to a right vector basis for the original C{t) by means of the 
transformation \v\i\t))) = S{t) \d\^^)) j which, applied to 
Eq. (IA7> . yields 

Lit) \V\^\t))) = \v\i-'\t))) + K{t) (A8) 

Equation ( IA8> exhibits an important feature of the set 
I (t))) |, namely, it implies that Jordan blocks are in- 
variant under the action of the Lindblad superoperator, i.e., 
the index a denoting the Jordan block is preserved under C. 

An analogous procedure can be employed to define the left 
basis. Denoting by j{{Qait)\ = {{Qa{t)\S{t) the left eigen- 
state of £j(i) associated to a Jordan block Ja, Eq. (IA2t leads 
to the normalized left vector 



((Qal 1^ = (0, ...,0,1) 



(A9) 



The additional left vectors are defined as (note that these are 
just the transpose of the right vectors in the Jordan basis) 



1,0,0, 



,0 



= 0,... ,0,1,0 



(AlO) 



which imply the following expression for the left basis vector 

{{£^\t)\ = j{{£^^\S-\t)fo,C{t): 

{{E^:\t)\ C{t) ^ {{£li+'Ht)\ + ((f « Wl A„W, (All) 
or, equivalently, 

C{t)^\E^:^))=\*a{t)\E^^)) + \El^+'^)). (A12) 

Here we have used the notation j{{£a^°' ^■'| = j{{Qa \ and 
j((£^"°^|^0. 

We can now derive the orthogonality and completeness re- 
lations. First, the left and right basis vectors are orthonormal: 

{iS^it)\Dfm = j{{E^^\S-\t)S{t)\v'-^^))j^5ap6^'. 

(A13) 

Second, it is clear that (since it is a standard ba- 
sis) the Jordan basis is complete in the sense that 
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Ea,/3;i,, [Df))j J{{£^^ I = I- Applying S[t) to the left and 
S~^{t) to the right of this equation we therefore find the com- 
pleteness relation 

\V^^\t))){{£i^Ht)\=I- (A14) 
As a final point of clarification, note that, even though in the 



stationary Jordan basis left and right basis states coincide (up 
to transposition), this is not the case in the time-dependent 
basis. This difference between left and right vectors is due to 
the non-unitarity of the similarity matrix S. To see this, note 
that for a 1 -dimensional Jordan block: \'Da))j = \£a)).j, so 
instead of {{£a{t)\ = j{{£a\S^^{t) we can write {{£a{t)\ — 
j{{Va\S-'{t). Then = {S-'Y it)\Va)) j, which 

does not equal \Va{t))) = S{t)\'Da)),j since S is not unitary. 
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